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1. INTRODUCTION 



The conventional heat pipe with its ultra high conductance has made 
great advances in solving many of the problems of temperature control. 
Basically, the heat pipe is a passive device - its operating temperature 
is determined by the source and sink temperatures and the heat load im- 
posed on it. Source temperature and heat load are coupled, i.e. , if the 
heat load varies, the source temperature will also vary. In many situa- 
tions, such as the cooling of sensitive electronic equipment, it is 
desirable to have a stable source temperature under conditions of varying 
heat load or sink temperature. The introduction of the variable conduct- 
ance heat pipe has provided a means of effectively achieving source temp- 
erature stabilization under conditions of varying heat load or sink temp- 
erature. Several techniques have been proposed for achieving variable 
conductance, but of all the proposed methods, the one which has received 
the most attention is the use of a noncondensible gas. 

Marcus [1]* has provided an excellent discussion of all the proposed 
control mechanisms and particularly the use of noncondensible gases. 

In recent years a great deal of work has been done in the development of 
gas loaded variable conductance heat pipes [2-8]. Most of this work, 
however, has been concerned with heat pipes operating in the zero gravity 
environment of space. During operation of a gas-controlled heat pipe, the 
flow of working fluid vapor tends to sweep the noncondensible gas toward 
the condenser end of the heat pipe where it collects and forms a gas plug. 
This plug of non-condensible gas blocks a portion of the condenser and 
very little vapor reaches it. Consequently, this portion of the condenser 
is unavailable for condensation and the effective heat transfer area of 



* Numbers in brackets indicate references listed in the Bibliography. 
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the condenser is reduced. If the thermal load to the heat pipe is in- 
creased, the effect is to increase the pressure and the mass-flow rate 
of the working- fluid vapor. This, in turn, compresses the noncondensible 
gas reducing the volume of the condenser which it occupies. This then 
increases the effective condenser area, thus accommodating the increased 
heat load. On the other hand, if the load on the pipe is reduced, the 
noncondensible gas will expand, thus reducing the effective condenser area. 

Inherent in the design of gas-controlled variable conductance heat 
pipes is the necessity for understanding the interaction between the work- 
ing-fluid vapor and the noncondensible gas. An adequate design implies 
the ability to predict the location of the interface, between the vapor and 
gas as a function of the heat ]oad, source, or sink temperature and any 
of the other parameters of the pipe which might be varied. 

The first analytical model which was used in the design of gas-con- 
trolled [9] heat pipes assumes that the interface between the vapor and 
the gas is very sharp. Axial conduction in the pipe wall and mass dif- 
fusion between the vapor and noncondensible gas are both neglected. Al- 
though the qualitative predictions of this so-called flat front theory 
are substantially correct, the quantitative results differ significantly 
from experimental observations. Specifically, the flat front theory 
underestimates the operating temperature of the heat pipe. The major 
source of error appears to be in neglecting axial conduction in the heat 
pipe wall. 

Rohani and Tien [10] have analyzed the two dimensional heat and mass 
transfer process in the gas-loaded heat pipe. Their analysis neglects 
axial conduction in the pipe wall but does include the effects of radial 
flow in the vapor space. 
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Edwards and Marcus [11] have formulated an improved analytical model 
to predict the performance of gas-loaded heat pipes. This analysis in- 
cludes the effects of axial conduction in the pipe wall, as well as one 
dimensional binary diffusion between the working fluid vapor and the non- 
condensible gas. The predications of this diffuse front theory seem to 
agree quite well with the experimental results obtained from a gas-loaded 
heat pipe using water as the working- fluid and air as the noncondensible 
control gas. 

The design procedures and analytical tools which have so far been 
developed to predict the performance of gas loaded heat pipes neglect the 
effects of gravity on the heat and mass transfer processes within the pipe. 
This leads to the tacit assumption that the vapor-gas diffusion process is 
for the most part one dimensional and that, except for the effects on 
capillary pumping in the wick, the diffusion process is independent of 
orientation. The experiments which are usually performed to confirm the 
diffuse front theory usually consist of measuring the axial temperature 
profile along the pipe surface by means of thermocouples. This procedure 
assumes that the surface temperature of the pipe gives a measure of the 
saturation temperature of the working fluid and thereby an indication of 
its partial pressure and concentration. The placement of the thermocouples 
along a single axial line assumes that the vapor-gas concentration profile 
within the pipe is one dimensional. 

The diffuse front theory has produced an extremely powerful tool for 
use in the design of gas-loaded variable conductance heat pipes. However, 
there still exist situations where one could be led to question the 
accuracy of the predictions of the diffuse front theory. The situation 
in which the molecular weights of the working fluid and the non-condensible 
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gas are very different is one in which effects of gravity must be given 
careful consideration. 

1.1 Objective 

The objective of this work has been to investigate experimentally 
the effects of gravity on the operation of gas-loaded, variable conduc- 
tance heat pipes in which the non-condensible gas and the working fluid 
are of significantly different molecular weights. 
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2. EXPERIMENTAL APPARATUS 



Three heat pipes have been built and tested. The dimensions of each 
heat pipe are given in table I. Each heat pipe is equipped with two 
pressure transducers , one at the evaporator end and one at the condenser 
end. The heat pipes are equipped with a precision bellows seal vacuum 
valve at the evaporator end. This is to facilitate loading with working 
fluid or noncondensible gas. At the condenser end each heat pipe is fitted 
with an "O-ring" sealed bleed screw. This is to facilitate the elimination 
of any extraneous noncondensible gases which might be brought into the 
heat pipe while filling with working fluid. Each heat pipe was mounted on 
a frame which would permit the pipe to be rotated through an angle of 180° 
so that the position of the pipe could be set from horizontal to vertical 
or any position in between. Each heat pipe has 24 surface thermocouples 
spaced at 5 cm. intervals along the adiabatic and condenser- reservoir 
sections. In addition, each heat pipe has two stainless steel sheathed 
thermocouples which penetrate 1.25 cm. into the evaporator and condenser 
of each pipe. The 2.5 cm. and 5 cm. diameter heat pipes have 6 extra 
surface thermocouples placed at 90° angular increments around the circum- 
ference of the pipe at two axial positions, 30.5 cm. and 76.2 cm. from the 
end of the condenser respectively. The heater assembly for each heat pipe 
consists of 3 mm. Nichrome resistance heater ribbon which has been painted 
with high temperature electrically insulating varnish and then wrapped 
around the full 30.5 cm. length of the evaporator section. The evaporator 
and adiabatic sections were insulated with several layers of thermal in- 
sulation. Thermocouples were placed between the layers of insulation to 
measure the heat loss through the insulation. Power to the heater is 
furnished by a regulated D.C. power supply. The current in the heater is 
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obtained by measuring the voltage across a stable precision resistor 
placed in series with the heater. All the instrumentation is read by 
means of a Hewlett-Packard 2010C Data Acquisition System. Further details 
of the construction of the experimental heat pipes can be found in the 
Theses of Humphries [12], Nayden [13], Batts [14], and Owendoff [15]. 

Figures 1,2 and 3 are photographs of the three heat pipes. 

2.1 Liquid Crystals 

The usual technique for experimentally investigating the transport 
processes taking place in the condenser-reservoir section of the heat pipe 
is to measure the axial surface temperature distribution along the outside 
of the pipe. It is assumed that the thermal resistance through the wick 
and the pipe wall is small enough that the surface temperature gives an 
accurate representation of the temperature of the working fluid vapor and 
noncondensible gas combination. It is also assumed that the working fluid 
is always at saturation conditions. The surface temperature is therefore 
assumed to give the local axial vapor saturation temperature and hence the 
saturation pressure. This saturation pressure corresponds to the partial 
pressure of working fluid vapor and is therefore representative of the 
local axial concentration of working fluid vapor. The use of thermocouples 
placed on a single axial line along the heat pipe tacitly assumes that the 
temperature and hence the concentration distributions are nearly one dimen- 
sional and vary only with the axial coordinate. 

To be able to measure any deviation from one dimensional behavior the 
outside surface of both pipes was coated with a mixture of microencapsu- 
lated cholesteric liquid crystals. Cholesteric liquid crystals are formed 
from esthers of cholesterol. These liquid crystals will exhibit all the 
colors of the visible spectrum as they are heated through their "event 
temperature range". Both the location on the temperature scale of the 
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event temperature and the width of the temperature range of the liquid 
crystals can be controlled by the proper formulation. An excellent com- 
prehensive review of the applications of liquid crystals to heat transfer 
and temperature measurement problems can be found in the recent work of 
Cooper, et al.[16]. The useful life of raw liquid crystals has been 
greatly extended by a process of microencapsulation. In this process, the 
crystals are coated with gelatin in a polyvinyl alcohol binder. This re- 
sults in the formation of small spheres of the order of 30 microns in 
diameter. In this form the liquid crystals are much better able to with- 
stand the deleterious effects of ultraviolet light and other environmental 
factors. The liquid crystals used in these experiments were in the form 
of a water based slurry. Three separate liquid crystals designated R-41, 
R-49, and R-56 were used. The R in the designation indicates that the 
width of the event temperature range is "regular”, i.e. on the order of 3°C 
from the onset of red to the fading of blue. The number in the designation 
indicates the temperature at which the colors begin to appear, e.g. R-49 
begins to show at red at about 49°C. Since the liquid crystals are viewed 
by reflected light, it is necessary to provide a black background against 
which to view them. 

Before applying the liquid crystals , both heat pipes were painted 
with a flat black enamel paint. Equal parts of the three liquid crystals 
(R-41, R-49, and R-56) were then mixed together. Several coats of this 
mixture were then applied to the 2.5 cm and 5.0 cm diameter heat pipes. 

Each coat was allowed to dry before applying the next. The 1.6 cm dia- 
meter heat pipe was coated with only the R-49 liquid crystal. During 
operation, the liquid crystal isotherms on the heat pipes were recorded 
by measuring the angular position at axial increments of at least 5 cm. 
of the green color line for a particular isotherm. 
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3 . EXPERIMENTAL PROCEDURE 



The filling procedure for each heat pipe consisted of first evacuating 
the heat pipe to a pressure of 10 microns of mercury by means of a vacuum 
pump. The proper amount of methanol or Freon was then added from a grad- 
uated bottle by regulating the flow with the bellows seal vacuum valve. 

This process always introduced a small amount of extraneous air. To 
eliminate the air sufficient power was applied to the heat pipe to raise 
the pressure inside to well above atmospheric. The heat pipe was allowed 
to come to a steady state and insure that all the air had accumulated at 
the condenser end. The bleed screw at the condenser end was then opened 5 
and the unwanted air was allowed to escape. The bleed screw was then 
closed. This procedure was repeated until the internal pressure at the 
condenser end was equal to the saturation pressure of methanol correspond- 
ing to the measured internal temperature at the condenser end. The par- 
tial pressure of air was then zero indicating all the air has been elim- 
inated. The heat pipe was then run as a conventional heat pipe to obtain 
its operating characteristics. The heat pipe was then filled with the 
noncondensible gas. The procedure for filling the 5 cm. diameter heat pipe 
with krypton will be outlined; the procedure is similar for the other gases 
and the other heat pipes. Initially with the pipe operating at 100 watts 
nominal power input, sufficient krypton was added to increase the pressure 
in the pipe by 93. kPa (13.5 psi). After the desired data was recorded, 
the pipe was vented through the bleed screw to remove approximately one 
third of the krypton by reducing the system pressure approximately 31. kPa 
(4.5 psi). The same procedure was followed to produce the third krypton 
load. The amount of gas in the pipe for a particular load was calculated 
from the pressure and temperature inside the pipe when the pipe was not 
operating and therefore at ambient temperature. 
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3.1 



Operating Procedure 



Each heat pipe was operated with at least three different amounts of 
each noncondensible gas. The gas loads for each heat pipe are summarized 
in table II. For each gas load heat pipe No. 1 was operated at power 
levels of from 10 to 50 watts in 10 watt increments. Heat pipe No. 2 was 
operated at power levels of from 20 to 50 watts in 15 watt increments. 
Heat pipe No. 3 was operated at power levels of from 25 to 150 watts in 
25 watt increments. At each power setting steady state data was recorded 
in both the horizontal and vertical positions. The power settings were 
subsequently corrected by subtracting the losses through the insulation 
to obtain the actual power to the heat pipe. 
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4. EXPERIMENTAL RESULTS 



The results for the largest and smallest diameter heat pipes exhibit 
the extremes in the effects of gravity on heat pipe operation: from 

essentially no effect on heat pipe No. 1 (1.6 cm. diameter) to very pro- 
nounced effects on heat pipe No. 3 (5 cm. diameter). The results for these 
two cases are therefore presented first followed by the results for heat 
pipe No. 2 (2.5 cm. diameter). 

4. 1 Heat Pipe No. 1 

Figures 4 and 5 show the measured surface temperature distribution 
from the thermocouples on heat pipe No. 1 (1.6 cm. diameter). They are plotted 
as the measured surface temperature minus the ambient temperature as a 
function of distance from the beginning of the condenser. These figures 
show the characteristic temperature distribution of the gas loaded heat 
pipe. The higher temperature region is indicative of the active condenser, 
the lower, near ambient temperature region is indicative of the gas-blocked 
inactive condenser and the intermediate region of changing temperature is 
indicative of the gas-vapor interface diffusion region. 

From Figure 4, with methanol and helium, the temperature profiles in 
the horizontal and vertical case are of the same shape , but the position 
has been shifted. The interface region is shifted toward the condenser in 
the vertical case, as compared to the horizontal. The lighter helium is 
tending to rise toward the higher condenser. With methanol and krypton, 
the shape of the profiles is altered somewhat. From Figure 5, a comparison 
of the horizontal and vertical profiles indicates that the vapor-gas inter- 
face region has been broadened in the vertical case, compared to the same 
power level in the horizontal case. This effect was more pronounced at 
the lower power levels. When operating in the vertical position the 
gravitational force causes the heavier krypton to tend to settle toward 
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the lower evaporator. This tendency is counteracted by the momentum flux 
of methanol vapor toward the condenser. The vertical operating position 
with methanol and krypton is a gravitationally unstable situation (heavier 
krypton above lighter methanol). This instability is borne out by the 
observation of the liquid crystal isotherm. In the vertical position the 
isotherm was in constant motion, randomly oscillating about a mean position 
with a period of the order of several seconds. The liquid crystal isotherm 
on heat pipe No. 1 was always nearly perpendicular to the axis of the heat 
pipe. This seems to indicate that the heat and mass transport processes 
within a heat pipe of this smaller diameter are nearly one-dimensional. 

4. 2 Heat Pipe No. 3 

Figures 6 and 7 show the surface temperature profiles on heat pipe 
No. 3. These results at first seem anomalous. In Figure 6, with methanol 
and helium the vertical temperature profile is consistant with what would 
be expected for a gas loaded heat pipe. The horizontal profile on the 
other hand is not. In the horizontal case, the pipe appears to have a 
temperature distribution which is continuously decreasing, indicating no 
clear regions of active and inactive condenser. In Figure 7 with methanol 
and krypton, both the horizontal and vertical cases have this indistinct 
temperature profile with no clear indication of a vapor-gas interface 
region. In this case, the liquid crystal isotherms are invaluable in 
interpreting the temperature distributions. 

Figures 8 through 19 are plots of the measured liquid crystal isotherms on 
heat pipe No. 3. The figures shown are plane projections of one side of 
the heat pipe. Results for all krypton loads and all helium loads for both 
the horizontal and vertical positions are shown. Figures 8, 9 and 10 are for 
two different amounts of krypton with the heat pipe in the horizontal 



19 



position. From these figures, it is obvious that stratification of the 
vapor and noncondensible gas has occurred. The heavier krypton has sunk 
to the lower portion of the heat pipe and the methanol vapor has risen to 
occupy the upper portion of the pipe. The isotherms indicate that the 
krypton is distributed along the full length of the condenser and even 
seems to penetrate into the adiabatic section.. This is particularly 
apparent in Figure 10 for the higher gas load. For both gas loads, as the 
power is increased, the stratification becomes more pronounced so that in 
Figure 9, at 94.1 watts the 56°C isotherm is nearly horizontal. 

Figures 11, 12 and 13 are similar plots for two different amounts of 
helium with the pipe operating in the horizontal position. In this case, 
the helium migrates to the top of the pipe and the methanol vapor settles 
to the bottom. Although less severe than with krypton, there is still a 
tendency toward horizontal stratification. As can be seen in Figure 13, 
at higher helium loads and higher power levels this tendency is increased. 

In all the horizontal cases, the liquid crystal isotherms indicate that 
the vapor transport and vapor-gas diffusion processes are highly three 
dimensional. This explains the seemingly anomalous surface temperature 
distributions obtained with the thermocouples. 

Figures 14 through 19 show the liquid crystal isotherms for the vertical 
position for both krypton and helium. In the vertical position the trans- 
port processes within the heat pipe are more nearly one dimensional with 
some modification depending on the particular gas. In Figures 14, 15 and 
16 , it is apparent that the gravitational force acting on the heavier kryp- 
ton has a tendency to spread out the vapor gas interface region. Indeed 
for some power levels and krypton loads, it appears to cover the full 
length of the pipe. It should be pointed out that the isotherms for the 
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krypton loaded pipe in the vertical position were in constant motion 
indicative of the gravitationally unstable nature of this case. 

The isotherms for the helium loads in the vertical position are 
shown in Figures 17, 18 and 19. In this case, gravity is aiding the 
methanol vapor momentum flux in separating the gas and vapor. The iso- 
therms are very closely spaced*, only their location is shifted as the 
power or helium load is changed. The behavior in the helium/ vertical 
case is in agreement with the usual assumptions for gas-loaded heat pipes. 
4. 3 Heat Pipe No. 2 

Heat pipe No. 2 (2.5 cm. diameter) was run using methanol as the work- 
ing fluid with either krypton or helium as the noncondensible gas. Addi- 
tionally, this pipe was also run using FREON - 113 as the working fluid 
and helium as the non- condensible gas. The combination of Freon with a 
molecular weight of 187.4 and helium with a molecular weight of 4 gave the 
largest difference in molecular weights between working fluid and non- 
condensible gas. It should be pointed out that no change in wick struc- 
ture was undertaken to accommodate the use of Freon as a working fluid. 

Figures 20 and 21 show the measured temperature distributions on heat 
pipe No. 2. From Figure 20, with methanol and helium the differences 
between the horizontal and vertical orientation are readily apparent, 
although not as pronounced as in Figure 6 for heat pipe No. 3 (5 cm. dia- 
meter) . With methanol and krypton the temperature distributions are as 
shown in Figure 21. Here, there appears to be no difference between the 
horizontal and vertical orientation. Again, by examining the isotherms 
obtained from the liquid crystals, as shown in Figures 22 to 31, it is 
apparent that the temperature distribution obtained from thermocouples 
placed on a single line along the axis of the pipe can be somewhat 
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misleading. Although these liquid crystal isotherms do not indicate 
the severe degree of stratification in the horizontal position that 
occurs with heat pipe No. 3 (5 cm. diameter), they do show a definite 
distortion of the temperature distribution. Prom Figures 22 to 25, 
for methanol and krypton, the effect of gravity seems to be to distort 
the isotherms because of the tendency of the working fluid and non- 
condensible gas to stratify horizontally. For methanol and helium, 

Figures 26 to 31, the tendency to horizontal stratification is less 
pronounced. The difference between the horizontal and vertical orien- 
tation seem to be a difference in the width of the diffusion interface 
region between the working fluid and the noncondensible gas. This is 
apparent from both the thermocouple data (Fig. 20) and the liquid 
crystal data. The results for heat pipe No. 2, as might be expected, 
lie between those for the larger and smaller heat pipes. 

A few experiments were conducted with the 2.5 cm. diameter heat 
pipe using Freon-113 as the working fluid and helium as the non-condensible 
gas. Figure 32 shows the axial temperature distributions obtained with 
the thermocouples. Here, also the differences between the horizontal 
and vertical positions are quite obvious. The liquid crystal data shown 
in Figures 33 to 36 exhibit the same qualitative features as were found 
in the methanol-helium data. In the horizontal position, with Freon as 
the working fluid, the liquid crystal isotherms exhibited a slight degree 
of assymetry. It was felt that this assymetry was probably due to non- 
uniformities in the wick structure which could have caused local varia- 
tions in the wick thermal resistance. 
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5. CONCLUSIONS 



It has been demonstrated that orientation in the gravitational field 
can have a serious effect on the operation of gas-loaded heat pipes in 
which the working fluid and noncondensible gas are of significantly dif- 
ferent molecular weights. For the smallest diameter pipe (1.6 cm.) the 
effect of gravity is to distort and displace the vapor-gas interface region 
while not seriously affecting the basically one-dimensional nature of the 
transport processes within the pipe. 

For the largest diameter pipe (5 cm.) the effects of gravity are more 
profound. In the horizontal operating position stratification of the 
working fluid vapor and noncondensible gas occurs. The vapor and gas 
separate into horizontal layers. In this case, the transport processes 
inside the pipe are highly three dimensional. In the vertical position 
with the large diameter pipe, gravity induces a distortion of the vapor- 
gas interface region which is more severe than for the smaller pipe. 

No attempt has been made to determine what role other factors such as 
capillary pumping or wall conduction might play in the observed phenomena. 
It has been assumed that the large molecular weight differences were the 
dominant factor. 
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Photograph of Heat Pipe No. 1 (1.6 cm Diameter) 






Photograph of Heat Pipe No. 
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(2.5 cm Diameter) 



FIGURE 2 
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Photograph of Heat Pipe No. 



3 (5 cm. Diameter) 



FIGURE 3 
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FIGURE 6. Surface minus Ambient Temperature vs. Condenser Length 

Methanol -Helium 
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Methanol -Krypton 
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Figure 8 . Liquid Crystal Isotherms - Horizontal - 2.75 gm Krypton 
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Figure 9. Liquid Crystal Isotherms - Horizontal - 5.26 gm. Krypton 
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Methanol with 7.39 gm Krypton 



Qnet 
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Methanol With 0.182 gm Helium 



Qnet = 71. 2W 
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Figure 12 . Liquid Crystal Isotherms - Horizontal - 0.29 qm. Helium 
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Figure 13. Liquid Crystal Isotherms - Horizontal - 0.63 gm. Helium 
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Figure Vi o Liquid Crystal Isotherms - Vertical 
2.75 gm. Krypton 
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Figure 15 Liquid Crystal Isotherms - Vertical 
5.26 9 m * Krypton 
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Figure 16. Liquid Crystal Isotherms - Vertical 

Methanol With 7.39 gm Krypton 
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Figure 17 # Liquid Crystal Isotherms - Vertical - 

Methanol With 0.182 gm Helium 
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Figure 18 . Liquid Crystal Isotherms - Vertical - 
0.29 gm. Helium 
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Figure 19. Liquid Crystal Isotherms - Vertical - 
0.63 gm. Helium 
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Qnet = 42. 5W 
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Liquid Crystal Isotherms - Horizontal - Methanol with 6.17 x 10 kg Krypton 



Qnet = 46. 1W 




48 



Qnet = 



18 . 2W 



32. 5W 



47. 3W 



Tamb = 24°C 



25°C 



27°C 



100 _ 



75 



jC 

-H 

b0 

C 

3 

d) 

w 

c 

d> 

TD 

C 

O 

O 



50 



25 





40^ 






45°_C^ 






53°C 

X 






40°C 

/ 






45°C 

iT 






^.53°C 





40°C 

45°C 

53°C 




Liquid Crystal Isothenns u - Vertical 
Methanol with 6 . 17 x 10 kg Krypton 

Figure 24 
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Liquid Crystal Isotherms^- Vertical 
Methanol with 1.59 x 10~ kg Krypton 

Figure 25 
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Liquid Crystal Isotherms - Horizontal - Methanol with 6.17 x 10 kg Helium 
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Liquid Crystal Isotherms,.- Vertical 
Methanol with 3.05 x 10 kg Helium 

Figure 29 
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Liquid Crystal Isotherms Vertical 
Methanol with 6.17 x 10 kg Helium 

Figure 30 
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Liquid Crystal Isotherms^- Vertical 
Methanol with 1.24 x 10 kg Helium 

Figure 31 
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Figure 32. Surface minus ambient temperature vs. 

Condenser Length; Freon - Helium 
(Heat Pipe No. 2) 
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Liquid Crystal Isotherms - Horizontal - Freon with 8.26 x 10 kg Helium 

Figure 33 
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Figure 35 
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